Detection of synchronised chaos in coupled map networks using symbohc dynamics 
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We present a method based on symbolic dynamics for the detection of synchronization in networks 
of coupled maps and distinguishing between chaotic and random iterations. The symbolic dynamics 
are defined using special partitions of the phase space which prevent the occurrence of certain 
symbol sequences related to the characteristics of the dynamics. Synchrony in a large network can 
be detected using measurements from only a single node by checking for the presence of forbidden 
symbol transitions. The method utilises a relatively short time series of measurements and hence is 
computationally very fast. Furthermore, it is robust against parameter uncertainties, is independent 
of the network size, and does not require knowledge of the connection structure. 
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Symbolic dynamics is a fundamental tool available to 
describe complicated time evolution of a chaotic dynami- 
cal system, the Smale horseshoe being a most famous pro- 
totype 0. In symbolic dynamics, instead of representing 
a trajectory by infinite sequences of numbers one watches 
the alternation of symbols. In doing so one 'loses' some 
information but certain important invariants and robust 
properties of the dynamics may be kept pj. Most stud- 
ies of symbolic dynamics in the literature are based on 
the so-called generating partition Q of the phase space, 
for which topological entropy of the system achieves its 
supremum j^]. On the other hand, finding a generating 
partition for general systems is a challenging open prob- 
lem and is the issue of the constant research j^. The 
misplacement of partition leads to 'diminishing' of the 
computed entropy of the system . Symbolic dynam- 
ics based on generating partitions plays a crucial role 
in understanding many different properties of dynami- 
cal systems. However, as we shall show, dynamics based 
on certain other partitions also have important practical 
uses. In the following, we construct specific partitions 
such that certain symbolic transitions do not appear in 
the time evolution of a dynamical system. We present 
two applications, namely distinguishing between chaotic 
and random iterations and the detection of synchroniza- 
tion in networks. In the first application, the absence 
of these forbidden transitions in a chaotic iteration dis- 
tinguishes it from a corresponding random iteration. In 
the second application, we utilise the fact that in a syn- 
chronized network all nodes follow the same evolution as 
an isolated unit, so that synchrony can be detected by 
comparing the symbolic dynamics from any single node 
against the dynamics in isolation. 

Synchronisation in large ensembles of coupled dynam- 
ical units is studied in many different fields fH 15 Ifll Hfl 
im H^ . There are several systems where synchronisa- 
tion processes are significant, e.g., information process- 
ing in the brain j^. The detection of synchronization 
of an extended system from local measurements has im- 
portant applications. For instance, certain pathologies 



in the neural system, such as epileptic seizures, manifest 
themselves by synchronized brain signals |lfij | and there 
is some evidence that they can be predicted by chang- 
ing levels of synchronization JTj. It is thus of interest 
to be able to determine different levels of synchroniza- 
tion in a brain area from local measurements, such as 
EEG recordings. Similarly, fast and easy detection of 
synchronization in communication or power networks is 
important for efficient management of such large systems 
|l2| . The methods so far known to detect synchronisation 
are based on the decrease in the topological entropy or 
on the correlation between the bivariate or multivariate 
time series observed at different nodes using tools from 
nonlinear time series analysis [isl . The first method 
requires the construction of phase space from time se- 
ries using delay coordinates 0, 0| , which itself is not 
a trivial task and typically requires extensive computa- 
tional effort Furthermore, its success depends on 
the availability of a sufhciently long time series, although 
it is not clear what length is sufficient 0|. Even with 
an infinitely long time series, the detection of synchro- 
nisation relies on the decrease of phase space dimension, 
and it is not clear whether any decrease is due to the 
synchronisation of nodes or to the unsynchronised low- 
dimensional chaos. Similarly, distinguishing chaotic and 
random systems is not always easy, and the issue is of 
both fundamental and practical importance 0|. In this 
Letter, we shall address these problems using symbolic 
dynamics based on specific partitions. 

We consider dynamical systems defined by the itera- 
tion rule 



(1) 



Where t e Z is the discrete time and / : — > S" is a 
map on a subset S of R". Let Si : i = 1, . . . , m be 
a partition of S, i.e., a collection of mutually disjoint 
subsets satisfying U-"jS'i — S. We assume that the Si 
are nonempty and m > 2 to prevent trivial cases. The 
symbolic dynamics corresponding to |^ is the sequence 
of symbols {. . . , St_i, St, St+i, • ■ • } where St = z if x(t) £ 
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FIG. 1: A partition showing avoiding sets. In the pairs 
(&, 53), (52, 54), (5*2, 5*5), (5*4, 55), (53, 5*1), (5*3, 52), (S's, Si), 
(85,82), {85,83), [85,84), (85,85), the first set avoids the 
second. 

Si. We say the set Si avoids Sj if 

n s, = 0. (2) 

Clearly, if Si avoids Sj, so does any of its subsets. We 
also talk about a self-avoiding set if (0) holds with i = j. 
The essence of our method is based on choosing a par- 
tition where Q holds for one or more sets in the parti- 
tion. The significance is that if Si avoids Sj, then the 
symbolic dynamics cannot contain the symbol sequence 
ij. The notion is extended in a straightforward way to 
the kth iterate of /. Thus, if n Sj = 0, then the 

symbol sequence for the dynamics cannot contain the 
symbols i and j at two positions which are k — 1 sym- 
bols apart. This constrains the symbol sequences that 
can be generated by a given map, and provides a robust 
method to distinguish between different systems by in- 
specting their symbolic dynamics. As examples of avoid- 
ing sets, we mention that for the common unimodal maps 
of the interval [0, 1], such as the tent or logistic maps, the 
set (x*, 1] and its subsets are self-avoiding, where x* de- 
notes the positive fixed point of /. More generally, if 
/ is a continuous one-dimensional function on the inter- 
val [ri,r2] having m fixed points xl, X2, --x^, and the 
derivatives / (x*) exist, then one can define a partition 
with Si = [ri,xl]. Si — {x*_j^,x*] for i = 2,...,m, and 
Sm+i — Then one can have several avoiding 

pairs, namely (Figure [T]l. 

if/'(a;*) > 0, then f{S^)nSj =9),yi < j 
if fix*) < 0, then f{S^) Sj = 0,Vi > j. 

Before turning to their applications, we first show that 
partitions which contain avoiding sets always exist and 
give a constructive proof. Suppose one starts with some 
partition of m sets for which does not hold for any 
i,j. Now fix some pair i ^ j- If f{Si) C Sj, then Si 
can be arbitrarily partitioned into two parts which avoid 
each other. If, on the other hand, f{Si) (f. Sj, then let 
Sm+\ = f~^iSj) n Si, and re-define Si as 5- = Si\Sjn+i 
to obtain a new partition with m + 1 sets, where S'i 
avoids Sj. The same argument can be used to construct 
self-avoiding sets: Assume f{Si) <f. Si (otherwise further 



partition Si to obtain a set which is not invariant under 
Z)!!! and define = r^{S,)r\S, and 5,' = S^S^^i, 
so that S[ is self-avoiding. Hence, it is possible to modify 
a given partition so that the sequence ij never occurs in 
the symbolic dynamics. 

We first consider the problem of distinguishing chaotic 
and random iterations. To be specific, suppose the iter- 
ation rule is given by the tent map 

/(x)^(l-2|x-l|). (3) 

Its stationary density on [0,1] is the uniform density, 
yix) = 1. If we take the generating partition, which is 
simply a partition of [0, 1] at the middle, then the above 
tent map and a random sequence from a uniform distri- 
bution give equal Kolmogorov-Sinai (KS) entropy [2^ . 
This can be seen by considering the permutation entropy 
based on transition probabilities of length d, 

PE = -- ^ p{ii,i2...id)Mp{h,i2,--;id) (4) 

ii,i2--- 

where p{ii, 12, ■ ■ • , id) is the joint probability that a;(l) G 
Si-^ , x{2) ^ Si^, . . . , x{d) € Si^ . Taking d = 2 and m — 2, 
the iteration takes a trajectory to the left or right parti- 
tion with equal probabilities, and the situation is identi- 
cal to the tossing of a fair coin. Hence, for this case the 
symbolic dynamics does not distinguish between the ran- 
dom and chaotic sequences, and the permutation entropy 
for both systems is equal to In 2. Now we shift the par- 
tition by S, i.e take the separation point at /p = 1/2 -t- S, 
and calculate the entropies of the tent map and the ran- 
dom iteration with the same probability distribution, as 
a function of 5. It turns out that the difference between 
the two entropies is maximum when Ip is chosen as the 
fixed point a/ (a+1) of the tent map. For the correspond- 
ing partition, i.c, S'i = [0,lp],S2 = (^p,l], the transition 
2 — > 2 does not occur. So, the random sequence can be 
identified by the presence of this transition. As the par- 
tition point is moved further to the right in the range 
a/{a + 1) < Ip < 1, the transition remains forbidden, 
but the difference between entropies decreases, reaching 
zero at /p = 1. The implication is that a longer time se- 
ries would now be needed as the observed occurrence of 
the symbol 2 becomes less frequent. Hence, the optimal 
partition would be the one for which the self-avoiding 
set ^2 is largest, which also yields the largest entropy 
difference between the chaotic and random maps. Note 
that the choice of Ip in the range (1/2, a/ (a -1-1)) gives a 
non-generating partition, but has no two-symbol forbid- 
den sequences. Although some longer sequences may be 
forbidden, their detection requires longer time series and 
more computational effort. On the other hand, avoiding 
sets, which forbid short sequences, is faster and can uti- 
lize shorter time-series. We next turn to the detection 
of synchronization in networks. Consider the well-known 
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FIG. 2: Detection of synchronisation in various networks of coupled maps. Figures (a)-(d) are plotted for coupled tent maps 
^ and (e)-(f) are for coupled logistic maps (42;(1 — a;)). The networks consist of (a) two coupled maps, (b) a globally coupled 
network of size 50, (c) a small-world network of size 100 and average degree 30, (d) a scale-free network of size 100 and average 
degree 15 constructed by the algorithm given in |2^. (e) a scale-free network of size 200 and average degree 10, and (f) a 
random network of size 200 and average degree 8. The horizontal axis is the coupling strength, and the vertical axis gives the 
synchronization measure (-) for the whole network, as well as the transition probability /3) (o) calculated using a scalar 
time series from a randomly selected node. In all cases, the synchronization region (cr^ = 0) coincides with the region where 
P{/3, beta) = and the other transition probabilities are nonzero, which is the situation for the uncoupled map. Note that in 
subfigure (c) there is an interval of e roughly between 0.65 and 0.8 such that P(/3, /?) — 0, but there is no synchronization as 
P{a, a) (shown by □) is also zero here, unlike the case for the isolated tent map. 



coupled map lattice model 

N 

x,{t+l) = /(a;,(<)) + ^^C., [f{x,{t))-f{xdm (5) 

where Xi{t) is the state of the ith node at time t, i — 
1, . . . , N, Cij are the elements of the adjacency matrix 
with value 1 or depending upon whether i and j are 
connected or not, e G [0, 1] is the couphng strength, and 
ki is the number of neighbours of the ith unit. Clearly, 
the system ((S)) can exhibit a much wider range of be- 
haviour than depending on the couphng structure and 
the choice of parameters; so the corresponding symbol se- 
quences observed from a node can vary widely. However, 
at the synchronized state Xi{t) = Xj{t) for all i,j and t, 
all nodes evolve according to the rule iQ). It follows that 
when the network is synchronized, the symbolic sequence 
measured from a node will be subject to the same con- 
straints as that generated by (P). If the partition used to 
define the symbolic dynamics contains avoiding sets, one 
has a simple way of detecting synchronization of the net- 
work by simply choosing a random node and observing 
the presence or absence of symbol subsequences which 
are forbidden in 

For numerical illustration, we take / to be the tent 
map with the partition Si = [0,a;*], 5*2 = (a;*,l], 
at the fixed point x* = a/ (a -I- 1), so that the transi- 
tion 2^2 has zero probability for the isolated map. 



In actual calculations, we use a slightly modified way to 
define the symbolic dynamics: To two consecutive mea- 
surements XtXt+i we assi gn t he symbol a ilxt+i > xt and 
the symbol (3 otherwise [l9|. This is completely equiva- 
lent to the symbol sequences defined by the sets Si , S2 
because of the specific partition point x* ; thus the tran- 
sition (3 ^ (3 does not occur for the single tent map. The 
advantage of using a, j3 is that one only needs to check 
increases and decreases in the measured signal, which 
adds more robustness in case the fixed point is not pre- 
cisely known. We evolve (|3J) starting from random initial 
conditions and estimate the transition probabilities using 
time series of length r = 1000 from a randomly selected 
node. Note that the length of the time series is much 
shorter time than would be required by standard time- 
series methods which use embedding to reconstruct the 
phase!^ space for large networks JJj. In our case, however, 
the length of the series is independent of the network size. 
We estimate the transition probability P{i,j) by the ratio 
X]t'T-(st = i,St+i = j)/ J2t''^(^t ~ where n is a count 
of the number of times of occurrence. Synchronisation is 
signalled when the variance of variables over the network 

given by = (^j^J^tl^ii^) ^''^P'^ *o ^^''O' 

where x{t) = jf'^iXi{t) denotes an average over the 
nodes of the network and (•■■)( denotes an average over 
time. Fig. 121 summarizes the results. It is seen in all cases 
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that the region for synchronization exactly coincides with 
the range for which (3) is zero (and other transi- 
tion probabiHties are nonzero; see subfigure (c)). Hence, 
regardless of network topology and size, both synchro- 
nized and unsynchronized behavior of the network can 
be accurately detected over the whole range of coupling 
strengths using only measurements from an arbitrarily 
selected node. 

We note that the method applies equally well to any 
diffusive network, say 

N 

x,{t + 1) = f{x,{t)) + w^jg{x,{t), xj{t)) (6) 
i=i 

which may be asymmetrically coupled or weighted, with 
weights given by Wij, x is n-dimensional, and g is a dif- 
fusion function satisfying g{x^x) = Va;. The basic idea 
is again the same: So long as the synchronized dynam- 
ics is identical to the isolated dynamics, synchrony can 
be detected by comparing the transition probabilities of 
symbolic sequences to those for the isolated map. The 
fact that certain transition probabilities are exactly zero 
makes this procedure especially robust, and the choice of 
the partition plays an important role here. For networks 
including transmission delays, the synchronized dynam- 
ics is no longer identical to that of the isolated map . 
Nevertheless, it can be shown that the properties of the 
map still induces certain forbidden transitions, whose 
presence or absence can be used to detect synchrony in 
the delayed network. Due to space constraints, we defer 
these results to an extended paper. 

In conclusion, we have presented a simple and effec- 
tive method based on symbolic dynamics for the detec- 
tion of synchronisation in diffusively coupled networks. 
The method works by taking measurements from as few 
as one single node, and can utilize rather short sequences 
of measurements, and hence is computationally very fast. 
Furthermore, it is robust against noise and parameter un- 
certainties, is independent of network size, and does not 
require knowledge of the connection structure. Addition- 
ally, the method allows to distinguish between chaotic 
and random iterations with same underlying probability 
density. 
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